Random-field Ising model on complete graphs and trees 
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We present exact results for the critical behavior of the RFIM on complete graphs and trees, both 
at equilibrium and away from equilibrium, i.e., models for hysteresis and Barkhausen noise. We 
show that for stretched exponential and power law distributions of random fields the behavior on 
complete graphs is non-universal, while the behavior on Cayley trees is universal even in the limit 
of large co-ordination. 



The central issue in the equilibrium random field Ising 
model (RFIM) is the nature of the phase transition from 
the ferromagnetic state at weak disorder to the frozen 
paramagnetic state at high disorder. The existence and 
universality class of the RFIM transition, is key as the 
best experimental tests of RFIM theory are diluted an- 
tiferromagnets in a field, which are believed to be in the 
same universality class as the RFIM [Q. After some con- 
troversy it was rigorously demonstrated that the RFIM 
transition occurs at a finite width of the distribution in 
three dimensions Q and at an infinitesimal width in one 
and two dimensions. Moreover, Aharony || showed that 
within mean field theory at low temperatures, the tran- 
sition is first order for bimodal disorder distributions but 
second order for unimodal distributions. Numerical stud- 
ies at zero temperature suggest that in four dimensions 
the bimodal case is first order and the Gaussian case is 
second order. The analysis in three dimensions is less 
conclusive ||. The difference between the Gaussian and 
bimodal cases has been attributed to percolative effects 
H . We have recently shown that at zero temperature, 
the mean-field theory is non-universal J(| in the sense 
that the order parameter exponent may vary continu- 
ously with the disorder. Exact optimization calculations 
[Q, ||] in three dimensions have also suggested that the 
correlation length exponent, as deduced from finite size 
scaling, is non-universal ||. 

Motivated by the fact that the RFIM is non-universal 
within mean-field theory for the stretched exponential 
distribution, we have analyzed the the RFIM on complete 
graphs with disorder distribution, (Sh/\h\) x (0 < x < 1, 
\h\ < Sh). We find that this distribution is anomalous 
in the sense that this sort of disorder never destroys the 
spontaneously magnetized state, at least within mean- 
field theory. The behavior of the RFIM on complete 
graphs is thus quite varied and anomalous. To deter- 
mine whether this non-universality extends to other lat- 
tices, we have analyzed the zero temperature RFIM on a 
Bethe lattice for the stretched exponential and power law 
distributions of disorder. We prove that the Bethe lat- 
tice is universal, provided the transition is second order, 
even in the limit of large co-ordination. This is surpris- 
ing since in this limit the Bethe lattice usually approaches 
the mean-field limit. 



We also extend the results outlined above to the non- 
equilibrium case. Ground state calculations of hysteresis 
and Barkhausen noise in the RFIM have demonstrated 
that the spin avalanches are controlled by the equilib- 
rium RFIM critical point |juj [uj. It is thus not surpris- 
ing, and we confirm, that the magnetization jump in the 
hysteresis loop is non-universal for the stretched expo- 
nential disorder distribution. The integrated avalanche 
distribution also has a non-universal exponent due to the 
non-universality of the order parameter. But the "differ- 
ential" mean-field avalanche exponent is universal even 
in cases where the order parameter exponent is not. In 
contrast, as expected from the equilibrium results, the 
Bethe lattice exhibits universal non-equilibrium critical 
behavior. 

The Hamiltonian of the random-field Ising model is, 
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where the exchange is ferromagnetic (Jij > 0) and the 
fields hi are random and uncorrelated. In the non- 
equilibrium problem we sweep the applied uniform field, 
H, from — oo to oo and monitor the magnetization at a 
fixed Jij = J and for a fixed disorder configuration {hi}. 
This model has been proposed as a model for Barkhausen 
noise by Dahmen et al. jnj. The local effective field re- 
sponsible for a spin-flip is 
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The condition for a spin to flip is that h\^ > 0. The ran- 
dom fields are drawn from a specified distribution p(h). 
To test universality, we use the following distributions 
which are defined on the interval -Sh < h < Sh, 
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and 
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We have shown that pi, which is the low field expansion 
of a stretched exponential disorder distribution, leads to 
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non-universality in the ground state of the equilibrium 
mean-field RFIM Here we extend that result to the 
non-equilibrium case. We then show that the distribu- 
tion p 2 destroys the RFIM phase transition, in mean-field 
theory and on trees, for —1 < y < 0. 

First we discuss the behavior of the ground state of the 
zero-temperature, mean-field RFIM. The magnetization 
is given by 
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were h c (m) 
tization is 
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-Jm — H. The energy at a given magne- 
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Extremizing with respect to the order parameter, m, 
yields the ground-state mean-field equation, 
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The non-equilibrium critical points are found from the 
susceptibility x = dm/dH, which from (Q) is given by, 
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The avalanche distribution, d(s,t) that gives the proba- 
bility of finding an avalanche of size s at parameter value 
t, is found using a Poisson statistics argument [0, which 
yields, 



d(s,t) ~ s- T e- 1 s = 8- T g(a a t) 
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where g(x) is a scaling function and t = 1 — 2Jp(Jm e + 
H). Experimentally, it is more natural to make a his- 
togram of all avalanches up to the critical applied field at 
which the magnetization changes sign. This "integrated" 
distribution behaves as, 



D(s, Sh) = s 
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where r = \Sh — Sh c \. For a Gaussian distribution of dis- 
order, j3 — 1/2, (J = 1/2, t = 3/2. We have shown, how- 
ever, that in the ground state for the distribution (||), the 
equilibrium order parameter exponent, j3 — 1/y. In con- 
trast it is evident from Eq. @ that the exponents a and 
t are universal. The non-universality in non-equilibrium 
behavior arises in the magnetization jump and in the 
shape of the non-equilibrium phase boundary, as we now 
demonstrate. Consider the distribution (||). Integrating 
(0) yields the mean field equations, 
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for Jm, + H > 0, and 

m = V-^-{Jm + H) + i \Jm + H\ y+ 1 (12) 

for Jm + H < 0. Here we have defined, J = J/8h, 
H = H/Sh. Setting H = in either (|lj) or © yields 
the equilibrium magnetization 
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At the critical point, the magnetization scales with the 
magnetic field as m e (r = 0,H) ~ H 1 / 6 . From Eq. @ 
it is evident that S = y + 1. The susceptibility x — 
dm/dH diverges when the barrier between the two local 
magnetization minima of the ground state energy ceases 
to exist. From (0), we have 



X = 



(y + l)[l - (Jm + H)y] 



y-(y + l)J[l-(Jm + H)y] 
and the critical condition 

y=(y+l)J[l-(Jm neq + H~ c ) y ]. 
This equation has the simple solution, 
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Substituting ( |16| ) into (|l l]), we find that the non- 
equilibrium magnetization jump is positive and has the 
value 
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for H — > Hj . Substituting this into Q16J) , the critical 
field is found to be, 
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This negative critical field is expected when starting 
with the positive magnetized state. By symmetry, the 
negative magnetization solution is at —H c . The value 
of the magnetization at that point is —m neq . Note that 
\m neq \ is not the size of the magnetization jump in the 
hysteresis loop. The jump in magnetization in the hys- 
teresis loop is Smhyst — \m neq \ +m(\H c \), where m(\H c \) 
is found by solving Eq. ([TJ]). The critical exponent asso- 
ciated with the jump in magnetization is determined by 
the behavior of the distribution p(h) at small fields, so 
that the critical exponents found here apply to distribu- 
tions of the form p(h) — exp(—(\h\/H) v ). For y < 1 these 
are the stretched exponential distributions ubiquitous in 
glasses, while for y > 2 they are more concentrated near 
the origin. 
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Figure 1: Top: The phase diagram of the non-equilibrium 
RFIM MFT using the disorder distribution (|), and Bottom: 
The magnetization jump at the phase boundary. In these 
figures, we took the exchange constant J = 1, The dotted 
line is for y = 0.5 while the solid line is for y — 2. Note that 
at the equilibrium critical disorder, Sh c , the hysteresis loop 
disappears. 

Now we briefly consider the distribution p2(h) given 
in Eq. (^). For y > this distribution is bimodal and 
it is easy to confirm the conclusion of Aharony || that 
the transition is first order. However the cases —1 < 
y < are more interesting. In these cases the disorder is 
dominated by small random fields, as the distribution is 
singular at the origin. It is easy to carry out the mean- 
field calculation (fjj) with the result, 

{8h\ 1+1,y 
meq = \ 1 ) Sh> J (19) 

By comparing the energies of E(m = 0), E{m = 1) 
and E(m eq ) (using Eq. @), we find that for Sh < J, 
the ground state is fully magnetized, while for Sh > J 
the ground state has magnetization (|l9|). The interest- 
ing feature of the result ( |l9| ) is that there is no phase 
transition at finite Sh, and the system is always ordered. 
The disorder distribution (0) thus destroys the ground 
state phase transition, due to the large number of small 
random fields. 

Now we determine whether the non-universal results 
found above for the mean-field theory extend to the 
ground state of the RFIM on a Cayley tree. The co- 
ordination number of a tree is taken to be z, while the 
probability that a spin is up is P+ and the probability 



Figure 2: Top The phase diagram for the non-equilibrium 
RFIM on a Cayley tree with coordination number a = 3 using 
the distribution (H) and taking the exchange constant J = 1. 
The dotted line is for y = 0.5 and solid line is for y = 2. The 
initial, linear part, of the phase boundary is due to the finite 
cutoff of the distribution (H). There is a discontinuity in slope 
of H c (5h) at the equilibrium critical disorder Sh c . Bottom 
The magnetization jump for the RFIM on Cayley trees for 
z = 4 and the distribution ([]), with the exchange constant 
J = 1. The dotted line is for y — 0.5 while the solid line is 
for y = 2. In both cases we find the same critical exponent, 
for example f3 — 1/2. In contrast, the mean-field result is 

P = l/y. 

that a spin is down is P_. The probability that a spin 
is up at level I can be written in terms of the probabili- 
ties at the level which is one lower down in the tree, this 
yields H H 

P +W = E ( "I ) P U l l)P-- g {l l)o+(a, fl ) (20) 

where a + (a,g) is the probability that the local effective 
field is positive when g neighbors are up. If we know the 
distribution p(hi) we can compute a+(a,g). Analyzing 
the equilibrium behavior, we have, 

/>oo 

a e _?{a,g) = / p{h)dh (21) 

The equilibrium Cayley tree model has been extended 
to the non-equilibrium case by considering a growth prob- 
lem in which the spin above the currently considered level 
in the tree is pinned in the down position ||l3|, Q. This 
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models the growth of a domain. The formalism is the 
same as in Eq. (EQ), with the modification that 



'(a,g) = a e ^(z,g). 



(22) 



From this equality and the form (|20|) it is easy to derive 
all of the non-equilibrium results from the equilibrium re- 
sults found using Eqs. ( po| ) and (|2~l|). To find the hystere- 
sis curve on a Cayley tree, we just shift the equilibrium 
magnetization as a function of field: by H — > H — J when 
sweeping from large positive fields and; by H — > H + J 
when sweeping from large negative fields. The behavior 
is evident in previous numerical work, but does not seem 
to have been noticed before. 

By direct iteration of the recurrence relation ( po| ) we 
show that a stable steady state solution, = 1 — P*, 
exists. It is easy to solve equation ( pcf ) in the steady state 
limit, at least for small values of a. For a = 1, 2 Cay- 
ley trees have no ordered state for any finite Sh, for the 
disorder distribution (||). But for a — 3 a ferromagnetic 
state does exist for a range of disorder. As we see from 
Eq. (po|), the a = 3 case leads to a polynomial of order 3 
which can be simplified to, 



rn 



— [mf(l -36 + a) - 1 + 3a + 36] = 



(23) 



were m = 2(P* - 1/2), a = a"; 9 (3,0) and b = a+ 9 (3, 1). 
Eq. ( |23l ) has the following solutions: 



to = 0: 



and m = ± 
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These solutions apply for any disorder distribution. For 
the distribution pi(h), performing the integrals yields, 



Ay - \2{y + 1) J + 3(3^ +1 + 1) f +1 

3(1 - 3y)T +1 
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We can now expand the magnetization around the critical 
point, J c , J = J c — (■ We find, 



(-12(y + l) + 3 



y+1 



{l + 3 y+1 )J V )e 
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(26) 

Thus to ~ e 1 / 2 for any y, so that (3 = 1/2 is universal. 
Since the non-equilibrium behavior on trees is related to 
that of the equilibrium behavior in such a simple manner, 
this universality extends to the hysteresis and avalanche 
exponents. It is easy to confirm numerically that the be- 
havior extends to large values of the branch co-ordination 
number a. Moreover by doing an expansion of ( p0[ ) using 
P+ = 1/2+TO, it is possible to show analytically that only 
the first and third order terms in to exist, regardless of 
the value of y in the disorder distribution (0) . This con- 
firms that for this distribution, the behavior is universal 
for all coordination numbers. 



For the distribution P2Q1) and a = 3 we get from 
Eq. @ 
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Just like we have done before we can expand to around 
the critical point, ,/ = J c — e: 



1/2 



Thus for p2 (3 — 1/2 is a universal exponent, too. 

In summary, on complete graphs (i.e. in mean-field 
theory) the RFIM at T — is non-universal. In particu- 
lar, the stretched exponential disorder distribution leads 
to a non-universal order parameter exponent and non- 
universal integrated avalanche exponent. In addition, the 
power law distribution has a regime in which a predomi- 
nance of small random fields destroys the transition and 
the RFIM always has a finite magnetization. In contrast 
the Cayley tree does not show either of these behaviors. 
Even in the limit of large coordination it is universal, 
with the usual mean-field order-parameter exponent 1/2. 
We have carried out some preliminary numerical studies 
of the behavior in three dimensions (with short range in- 
teractions) and find that the power law distribution of 
random fields does not destroy the transition. Moreover, 
the exceedingly small value of /3 in three dimensions ren- 
ders any non-universality in (3 a moot point. However the 
behavior in dimensions higher than three, or for longer 
range interactions in three dimensions could be more in- 
teresting. Finally, even for short range interactions in 
three dimensions, there have been suggestions of non- 
universality in the finite size scaling behavior ||. It is 
unclear, as yet, whether that behavior is related to the 
non-universality seen here. 

This work has been supported by the DOE under con- 
tract DE-FG02-90ER45418. 



Electronic address: 



dobrin@pa.msu.edu 



[1] S. Fishman and A. Aharony, J. Phys. C 12, L729 (1979). 

[2] J. Z. Imbrie, Phys. Rev. Lett. 53, 1747 (1984). 

[3] A. Aharony, Phys. Rev. B 18, 3318 (1978). 

[4] M. R. Swift, A. J. Bray, A. Maritan, M. Cieplak, and 

J. Banavar, Europhys. Lett. 38, 273 (1997). 
[5] M. R. Swift, A. Maritan, M. Cieplak, and J. R. Banavar, 

J. Phys. A 27, 1525 (1994). 
[6] P. M. Duxbury and J. H. Meinke, Phys. Rev. E 6403, 

6112 (2001). 

[7] H. Rieger, Frustrated Systems: Ground State Properties 
via Combinatorial Optimization, vol. 501 of Lecture Notes 
in Physics (Springer Verlag, Heidelberg, 1998). 

[8] M. J. Alava, P. M. Duxbury, C. F. Moukarzel, and 
H. Rieger, Phase Transitions and Critical Phenonema 
Volume 18 (Academic Press, 2001). 



5 



[9] J. A. D'Auriac and N. Sourlas, Europhys. Lett. 39, 473 
(1997). 

[10] K. Dahmen and J. P. Sethna, Phys. Rev. B 53, 14872 
(1996). 

[11] J. P. Sethna, K. A. Dahmen, and C. R. Myers, nature 
410, 242 (2001). 



[12] R. Bruinsma, Phys. Rev. B 30, 289 (1984). 
[13] D. Dhar, P. Shukla, and J. P. Sethna, J. Phys. A 30, 
5259 (1997). 

[14] S. Sabhapandit, P. Shukla, and D. Dhar, J. Stat. Phys. 
98, 103 (2000). 



